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We present power-law relaxation behavior of the local magnetizations in the equilibration dynam¬ 
ics of the spin-1/2 Heisenberg spin chain as an isolated integrable quantum system. We perform 
the exact time evolution of the expectation values of the local spin operators by evaluating them 
with the determinant formula of the form factors. We construct such an initial quantum state that 
has a localized profile of the local magnetizations, and perform the exact time evolution over a 
very long period of time. We show that the local magnetization relaxes as some power of the time 
variable with no definite time scale, while the fidelity relaxes very fast with its relaxation time being 
proportional to the inverse of the energy width, i.e. the Boltzmann time. 
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INTRODUCTION 


The non-equilibrium dynamics of isolated integrable 
quantum systems such as quantum quench has attracted 


much attention, recently |l|-|15|. Due to novel exper¬ 
imental demonstrations of out-of-equilibrium dynamics 
of closed quantum systems made of cold atoms fl6l4l8l | . 
it has become an extensive research subject how ther- 
malization occurs in isolated quantum systems For 
integrable quantum systems, which have a large num¬ 
ber of conserved quantities, equilibration dynamics is not 
trivial. It has been conjectured that integrable quantum 
systems after long time the expectation values of local 
quantities approach those of the Generalized Gibbs En¬ 
semble (GGE) Q- It has been investigated for particular 
integrable models and initial states whether the GGE 
conjecture holds or not 0-E3- 


The studies of non-equilibrium quantum dynamics are 
closely related to the recent renewed interest in the foun¬ 
dation of quantum statistical mechanics One of 

the fundamental questions is whether and how equilib¬ 
rium distribution functions are realized only through the 
unitary time evolution of a quantum state in an isolated 
system [27J. It has been demonstrated that in some gen¬ 
eral or concrete settings a pure initial state of an isolated 
quantum system equilibrates or thermalizes in a certain 
mathematical sense EMI 0, EMU. It is related to the 
idea that a typical pure state of a macroscopic quantum 
system can fully describe thermal equilibrium 22H24J . I26| . 
The question is often discussed in terms of the eigenstate 
thermalization hypothesis [32 . 


Although the above mentioned studies show that a 
wide class of isolated quantum systems thermalize or 
equilibrates to GGE, the time scale for thermalization 
has been studied only very recently El El- Further¬ 


more, the quantum dynamics of a pair of spins have been 
per formed numerically f34L [35| in an experimental setting 
[l8l |. However, it has not been explicitly shown how the 
expectation value of a certain local operator of an inter¬ 
acting integrable quantum system such as the Heisenberg 
chain of finite size equilibrates in time. It has not been 
clear even whether it approaches some constant value or 
not. Here we do not take the time average of the ex¬ 
pectation value. It is natural to conjecture that as the 
system size increases, the expectation value of a local 
operator at a given time becomes close to the limiting 
value of the infinite system. However, we expect that 
there are many aspects in the way how it approaches the 
equilibrium value in time and they are related to interest¬ 
ing topics of non-equilibrium dynamics close to integrable 
systems such as prethermalization j3(| and pre-relaxation 


37], which should occur between the integrability and its 


breaking under non-integrable perturbations. Moreover, 
it is important to give a concrete example of equilibration 
of a local operator in a long period of time for interact¬ 
ing quantum systems from the viewpoint of stability and 
equilibrium properties of KMS states [38|. 

Let us introduce the Hamiltonian of the anti¬ 
ferromagnetic Heisenberg spin chain under the periodic 
boundary conditions (PBC): crff +1 = of for a = X, Y, Z. 


n = \Y, {°f°j+i + a J a J +1 + (°f°f+1 - !)} ■ (!) 

3 = 1 

Here, er“ denote the Pauli matrices. We also call the sys¬ 
tem the XXZ chain. In the present paper we show the 
exact relaxation dynamics of local magnetizations (c * % ) 
of the spin-1/2 Heisenberg chain for an initially localized 
quantum state. We also study it for other states con¬ 
structed from spinon states. Although it shows a strong 
oscillating behavior in time, we argue that the amplitudes 
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of oscillation show a power-law decay, by evaluating the 
square deviations of the local magnetizations. We sug¬ 
gest the power-law decay may be universal for the expec¬ 
tation values of other local operators while the decaying 
exponent depends on the initial state. 

The paper consists of the following. We briefly review 
the Bethe-ansatz equations and the scheme to evaluate 
local magnetizations, which is due to the recent devel¬ 
opment of the algebraic Bethe-ansatz [39l l4(i| . We first 
show that the relaxation time of the fidelity is given by 
the Boltmann time and is consistent with recent rigorous 
results. Here, the fidelity is given by the square ampli¬ 
tude of correlation between the initial state |SB(0)) and 
state at time t: F(t) = |('L(t)|'I , (0))| 2 . We then show 
the time evolution of the profile of local magnetizations 
and observe how initially localized profile collapses 
in time. It seems that the resulting oscillations of local 
magnetizations ( 0 ^) continue forever. However, we shall 
show that the square deviations of the local magnetiza¬ 
tion exhibits a power-law decay. 


METHOD 


The eigenvectors and eigenvalues with M down spins 
of the Hamiltonian are constructed from solutions of the 
Bethe-ansatz equations (BAE) [dll . |42[ . 

2tt 1 M 

2 tan -1 (2A a ) = ^-J a + V 2 tan -1 (A« - \p), 

iV /3=1 

for a = 1,2,..., M. (2) 


Here, quantum numbers J a are given by integers or half¬ 
integers according to the condition: J a = (N — M + l)/2 
(modi). Here we take the branch: |tan -1 x| < 7r/2 for 
any x £ R. We denote by |{Aj}) the eigenvector as¬ 
sorted with a solution of the BAE, Ai, A 2 ,..., Am- The 
energy eigenvalue E of the eigenstate | {Ay }) is given by 


Ex 


M 


E 


2 

1 + 4AT 


(3) 


We consider a state |<f>) consisting of a superposition 
of the Bethe-ansatz eigenvectors | { Ay }) 


I*) = E|A> c a*- ( 4 ) 

Aes 


Here we have abbreviated eigenstate |{Ay}) by |A) and E 
denotes a certain set of solutions of BAE. 

We denote by (cr„(t)) the local magnetization on the 
mth site at time f, and define it by the expectation value 
of the local spin operator < 7 ^ with respect to the quantum 
state |<f>) at time t as follows. 


«(*)) = <$(t)| <|$(t)). (5) 


fid. 



FIG. 1: Time evolution of fidelity F(t) for states $(0)) given 
by sums of spinon states of equal weight in energy width A E 
for AE = 0.01, 0.05 and 0.1, respectively, in the XXX chain 
of N = 1, 000 in the sector of M = 499 down-spins. 


Here the time evolution of the state |$(t)) is given by 
|$(i)) = exp(— iHt/h)\$>(0)) where |$(0)} is given by |d>). 
We evaluate the local magnetization by the form factor 
expansion: 

<*wi <\m) 

= E^l cr ^l A ^^^^ A l 4 'l >exp ( _ ^- £; A)i) . (6) 

\,n 


Here we recall that E\ denotes the energy of the Bethe 
state |A). 

In order to evaluate the matrix elements (/z|cr^JA) of 
the local operator cr^ between two Bethe eigenstates |/x) 
and |A), we make use of the determinant formula of the 
matrix elements of tr^ derived by Kitanine, Maillot and 
Terras [39j |. For numerical evaluation, we further divide 
it out by the Cauchy determinant similarly as that of the 
ID Bose gas [4p[ . 

We remark that in order to show equilibration of local 
observables systematically, we may consider the Bloch 
vector, which is given by the expectation values of the 
traceless operators in the Hilbert space of the XXX spin 
chain with respect to a given quantum state [a|. We 
consider the local magnetizations at sites m, (tr^) for 
m = 1, 2,..., TV are elements of the Bloch vector. 


FIDELITY 


Let us evaluate the fidelity for pure states |<I>(0)) which 
are given the sums over spinon states with equal weight 
in the energy width A E for A = 0.01, 0.05 and 0.1, etc., 
specifically. Each spinon eigenstate |A) corresponds to a 
real solution of BAE |2]) specified by putting two holes in 
the set of quantum integers Jj s [4l|, |42j. Here we assume 
that the system size is given by N = 1,000, and the 
number of down-spins M by half the system size minus 
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FIG. 2: Relaxation time t = Tp of fidelity for the states 
given by sums over spinon states with energy width A E and 
the line: Tp = h/AE. 


1 : M = 499. 


I*(0)) = 5>.|A> 

a es 


(7) 


where symbol S denotes a subset of the set of all spinon 
eigenstates, and in the case of equal weight we set = 
1/VJ5[, where |S| denotes the number of elements in the 
set S. We plot in Fig. |T|the time evolution of the fidelity 
for the three states |$(t)) after initial time t = 0 with 
A E for A = 0.01, 0.05 and 0.1: F(t) = |(<F(t)|$(0))| 2 . 

We observe that the three graphs of fidelity F(t) versus 
time t shown in Fig. [1] are well approximated by the 
following expression [431 ] 


m = 


f. 


2 h 


t 2 + (/3h) 2 \AEj Sm “ l 2 H t 


f 


A E 


( 8 ) 


There is only one fitting parameter A E, practically. Here 
we assume that the inverse temperature is given by the 
energy width: 1//3 = A E. The initial Gaussian behavior 
should be consistent with other approaches 44, 4f| . Non¬ 
exponential behavior in the fidelity of interacting many- 
body systems is recently argued [461 ]. 

We estimate the relaxation time of fidelity through the 
numerical plots in Fig. [Q In each graph we determine it 
by the time when the fidelity takes the minimum value 
first after initial time. The red dots in Fig. |T]show the 
points of time when we determine the relaxation time 

t f . 

The estimates of relaxation time of fidelity Tp are 
plotted in Fig. [2] They are consistent with the line: 
Tp ft: h/AE. Here we remark that in the units of 
the present paper we take h = 1. Thus, it is given 
by the Boltzmann time assuming kpT = AE. Fur¬ 
thermore, it is interesting to observe that it coincides 
with a rigorous theoretical estimate [29] of thermaliza- 
tion time in a macroscopic quantum system for a typical 
non-equilibrium subspace: Tp ~ h/AE. 

For the XXX chain of N = 10 we evaluated the fi¬ 
delity for considering all the Bethe-ansatz eigenstates as¬ 
sociated with not only real solutions but also complex 
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FIG. 3: Collapse of an initially localized magnetization pro¬ 
file {(Tm(t)) (m = 1,2,..., N) in time evolution for the XXX 
chain with N = 50 in the sector of M = 24 down-spins. The 
profile is localized around at m = 25 when t = 0. 


solutions as in Ref. [53] (see also [48j|f. We compared 
the time evolution of fidelity for the state consisting of 
only real solutions and that for the state of all solutions 
including complex solutions. However, there is no clear 
difference in the time evolution of fidelity between the 
states of only real solutions and those of all solutions. 


LOCAL MAGNETIZATION 

Let us now consider mainly the all-spinon state l'F(O)), 
which is given by the sum over all the spinon states for the 
XXX spin chain with an even system size N in the sector 
of M = TV/2 — 1 down-spins, i.e. in the case of an almost 
half-filled lattice. We evaluate the local magnetization 
(< 7/n(t )) at time t by eq. (0. 

We observe in Fig. [3] that the profile of local mag¬ 
netizations (cr^(0)) for m = 1,2,...,TV at initial time 
t = 0 is localized around at the middle site of m = N/2 
and almost vanishes on other sites. We may consider it 
as the “positive” analogue of a “quantum dark soliton” 
constructed in the one-dimensional Bose gas with delta- 
function interactions ]3l], which are given by the sum 
over the type-II excitations in a branch. Here we remark 
that the average of a ^ is given by (TV — 2M)/N = 2/TV = 
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FIG. 4: Slow relaxation of local magnetization at site 1, 
(erf (£)), for the sum of all-spinon state in the XXX chain of 
N = 50 with M = 24. 


0.04 for N = 50. 

There are several procedures at different time steps as 
shown in Fig. [3] for the whole time evolution in which the 
initially localized profile of local magnetizations (cr^(t)) 
collapses to a uniformly oscillating profile in the final 
stage. First, the initially-localized profile separates to 
two localized waves moving in different directions. Then, 
the two localized waves propagate and collide each other 
at t = 18 ~ 20 around at the 1st site. They further prop¬ 
agate with smaller height of localized waves, and finally 
they merge into such a profile that is roughly uniform in 
space but strongly oscillating in time. 

Thus, the initially localized profile in Fig. [3] dynami¬ 
cally collapses to the roughly uniform profile oscillating 
in time. Here, the exact time evolution of the local mag¬ 
netization at a site shows strong oscillating behavior in 
Fig. [4] It seems that the fluctuations of the local magne¬ 
tizations do not easily relax to the order of the inverse of 
the number of the Bethe-ansatz eigenstates consisting of 
the quantum state 141), i.e. 0(1/|5|), which is the order 
of magnitude expected statistically. 

However, by taking the site-average of the squared de¬ 
viations of the local magnetizations over all sites we ob¬ 
serve a power law decay over a long period of time as in 
Fig© We define the site-average of squared deviations 
of local magnetizations on the mth sites (u^) by 

i N ( i N 

<(A^(i)) 2 ) = Jj £ ««} - 

m=l Y 3=1 

(9) 

Here we remark that due to the symmetry of the Hamil¬ 
tonian the sum of all the local magnetization is given by 
a constant: J2m=i a m = N — 2 M. 

For the all-spinon state of N = 30 , the squared devia¬ 
tions of local magnetizations decay almost as an inverse 
of time initially (0.0289/t 1 ’ 02 ), and as an inverse power 
of time with a smaller exponent (0.0765/t 0,806 ). For the 




FIG. 5: Relaxation of the site-average of the squared devia¬ 
tions of the local magnetizations © for IV = 50. Power law 
behavior such as 0.0250/t°' 68 and 0.0176/i°' 89 appears. 


all-spinon state of N = 50 the squared deviations of local 
magnetizations are shown in Fig. [5] For various other 
states which are given by the sums over sets of spinons, 
we observed similar power-law decay behavior of squared 
deviations ©. Here, the estimates of the exponent are 
different among the states. For instance, the yrast state, 
which is given by the sum over lowest excitations among 
spinon states, the estimate of the exponent is close to 
zero. We suggest that power-law decay may be universal 
in the time evolution of the expectation values of local op¬ 
erators for interacting integrable quantum systems, while 
the exponent depends on the states. 

We now argue that the time scale for equilibration of 
the local magnetizations is very long. We suggest that the 
time scale of the collapse of the initially localized mag¬ 
netization profile in the first stage is given by the system 
size divided by the spinon velocity. It is proportional to 
the system size. Here we recall that two localized waves 
travel in time in opposite directions in Fig. © However, 
the power-law decay of the squared deviations of the lo¬ 
cal magnetizations shown in Fig. [5] shows that there is 
no definite relaxation time in the whole long-time equili¬ 
bration dynamics of the local magnetizations. 

In summary we have performed exact relaxation dy¬ 
namics of the local magnetization profile for the initially 
localized state l^) in the XXX chain. Through the 
square deviations of local magnetizations we presented 
the power-law decay of the local magnetization in a long 
time. We suggest that the expectation values of other 
local operators may also decay as a power of the time 
variable for interacting integrable quantum systems. 
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